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Abstract. Despite the significance of type la supernovae (SNela) in many fields in 
astrophysics, SNela lack a theoretical explanation. We investigate the potential contri- 
bution to the SNela rate from the most common progenitor channels using the binary 
population synthesis (BPS) code SeBa. Using SeBa, we aim constrain binary processes 
such as the common envelope phase and the efliciency of mass retention of white dwarf 
accretion. We find that the simulated rates are not sufficient to explain the observed 
rates. Further, we find that the mass retention efficiency of white dwarf accretion sig- 
nificantly influences the rates, but does not explain all the differences between simulated 
rates from different BPS codes. 



1. Introduction 

SNela are generally thought to be thermonuclear explosions of carbon/oxygen (CO) 
white dwarfs (WDs ). The standard scenarios involve white dwarfs reaching the Chan- 
drasekhar mass (e. g. INomotolll 982 ) : either by hydrogen accretion fro m a non-degenerate 



companion (single-degenerate channel, SD; Whelan & Iben 1973) or by a merger of 
two CO WDs (double-degenerate channel, DD; W ebbink. 1984 ; Iben & Tutukov.1984) . 

We investigate the contribution from the DD an d SD channel to the SNIa rate with 
the binary popu lation synthesis (BPS) codes SeBa (Portegies Zwart & Verbuntlll996l; 



Nelemans et al.ii2001.. and Toonen et al. submitted). BPS codes are very useful tools 



to study the evolution of binary stars and the processes that govern them. We study 
the SNIa delay time distribution (DTD), where the delay time is the time between the 
formation of the binary system and the SNIa event. In a simulation of a single burst 
of star formation the DTD gives the SNIa rate as a function of time after the starburst. 
The DTD is linked to the nuclear timescales of the progenitors and the binary evolution 
timescales up to the merger. 
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Figure 1 . Simulated distribution of population of visible double white dwarfs as 
a function of orbital period and mass ratio, where mass ratio is defined as the mass 
of the brighter white dwarf divided by that of the dimmer white dwarf. Left model 
■ya is used, on the right model aa. The intensity of the grey scale corresponds to 
the density of objects on a linear s cale. The sam e grey scale is used for both plots. 
Observed binary white dwarfs (Marsh etal.ll2011b are overplotted with filled circles. 



2. Double degenerate channel 

We set out to predict SNIa rates for the double degenerate channel, with the additional 
constraint that our model corresponds well to the observed population of close double 
WDs (cDWDs) - of all flavours and masses. Even though there are no certain DD SNIa 
progenitors among the observed cDWDs, the DD SNIa and observed cDWD progeni- 
tors have gone through similar evolution paths and are affected by the same binary and 
stellar processes. 

cDWDs are believed to encounter at least two phases of common envelope (CE) 
evolution. In spite of the importance of the CE phase, it remains poorly understood. 
Severa l prescriptions for CE evolution have been proposed. The a formali sm ( WebbiiiS 
( 1984 )) is based on the conservation of orbital energy and the y formalism (iNelemans et al. 



(120001) ) is based on the conservation of angular momentum. In model aa the a formal- 
ism is used to determine the outcome of every CE. For model ya the y-prescription 
is applied unless the binary contains a com pact object or the CE is triggered by the 



Darwin-Riemann instability (lDarwinlll879l) . The ya-model reproduces the mass ratio 



distribution of the observed double white dwarfs best, see Fig. [T] 

The DTD of model ya and model aa are similar in showing strong declines with 
time and comparable time-integrated rates of 2.0 • 10"^ Mq^ resp. 3.3 • lO"'^ ^q^- 
Most importantly, the simulate d time-integrated numbers do not match the observed 
number of 2.3 ± 0.6 • 10"^ M'^ dMaoz etal.ll201ll ) by a factor of ~ 7 - 12. Many things 



influence the normalisation of the SNIa rates; the assumed binary fraction, metallicity, 
initial distribution of masses and orbital parameters. However, preliminary results show 
that the integrated rates are not affected by factors sufficient to match the observed 
rate. If so, the main contribution to the SNIa rate comes from other channels as for 
example the single degenerate scenario (e.g. supersoft sources), double detonating sub- 
Chandras ekhar accretors (see e .g. Kromeretal. 2010) or Kozai oscillations in triple 



systems dShappee & Thompson! (120121) : Hamers et al. in prep.). For more information 



about our study, see Toonen, Nelemans, Portegies Zwart submitted. 
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Figure 2. SNIa rate in the DD channel per yr per Mq formed stellar mass as a 
function of delay time. Delay times are shown for two different prescriptions of the 
CE phase; model -ya and model aa, see section [2l Overplotted wi th black circles are 
the observed values of the SNIa rate (see lMaoz & Mannuccil20ni, for a review). For 
comparison the grey circles show the observations scaled down by a factor 10. 



3. Single degenerate channel 



Theoretical rates for the SD channel predicted by different BPS codes vary over four or^ 



ders of magnitude and do not agree with each other or with observational data ( Nele mans et al. 
2012h . The exact origin of these differences remains unclear. We therefore study the 
effect of the efficiency of WD accretion, which is poorly understood beca use of pro- 
cesse s such as novae and stable burning. Three prescriptio ns are used (Nomoto et al. 



2007: Ruiter et al 



l2009l: lYungelsonl |20 lOl) that are based on iNomoto et all (l2007h and 



Prialnik & Kovetzl (Il995 ^) in different combinations which differ strongly, see Fig. [3^. 
We use the BPS code SeBa to simulate the SNIa rates in the SD channel for each reten- 
tion efficiency. The simulated SNIa rates are significantly affected by the prescription 
used for the mass retention efficiency. The integrated rates vary between 7.0- 10~^ Mq^, 
2.2 ■ 10"^ M^^ and an upp e r limit of 1 - 1 0"^ M ^^ wh en the reten t ion effi ciency is used 



as m 



Nomoto et all (l2007h . iRuiter et all (|2009[) rest ). lYungelsonl (l2010h where the ob- 
served rate is 2.3 ± 0.6 • 10"^ (Maoz et aLll201lb . However, the retention efficiency 



does not explain all differences between the theoretical SNIa rate distributions from 
different codes. For more information about this study, see Bours, Toonen, Nelemans 
in prep. 



4. Outlook 

To understand the differences in the predictions of the various BPS codes in the SD and 
DD scenario, we started a collaboration to comp are three BjPS co des. The codes in- 



volved are the Binary _c cod e (Izzard et al. 2006; Claeys et al.ll201 li) . the Brussels code 



( De Ponder & Vanbeverenl; [Mennekens et al. 2010) and SeBa. The comparison focuses 



on ffie evolution of low and intermediate mass binaries containing one or more white 
dwarfs. The goal is to investigate whether differences in the simulated populations are 
due to numerical effects, or whether they can be explained by differences in the input 
physics. The comparison indicates that the differences are caused by varying assump- 
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Figure 3 . a) The mass retention efficiency of WD hydrogen accretion as a function 
of the mass transfer rate. ?7tot represents the fraction of the transferred mass that is 
retained by the WD. If //tot = no mass is retained by the WD, when //tot - 1 the 
WD accretes matter most effectively. In the figure we have assumed the white dwarf 
to be of one solar mass, b) The delay time distribution that result from the different 
retention efficiencies discussed in the text. Note that we found no SNIa when the 
retention efficiency of iYungelsoni (i201ft) are used. 

lions for the input physics. For more information about this study, see the contribution 
to this volume by Claeys et al. 
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